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We consider the motion of conductivity spectral weight at a finite-temperature phase transition 
at which d^2_y2 density-wave (DDW) order develops. We show that there is a shift of spectral 
weight to higher frequencies if the quasiparticle lifetime is assumed to be isotropic, but a shift to 
lower frequencies if the quasiparticle lifetime is assumed to be anisotropic. We suggest that this 
is consistent with recent experiments on the pseudogap phase of the cuprate superconductors and, 
therefore, conclude that the observation of a downward shift in the spectral weight at the pseudogap 
temperature does not militate against the DDW theory of the pseudogap. 
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I. INTRODUCTION 

It has been proposedi that the anomalous behaviors 
observed in the pseudogap phase of the underdoped 
cuprates can be explained by assuming the existence of a 
new long-range ordered density- wave with dx^-y^ sym- 
metry (DDW)^ which competes with superconducting 
order. As a result of this order, the excitation spec- 
trum acquires a partial gap with the same symmetry 
and structure as the gap observed in the pseudogap 
state. The DDW theory has been applied to explain sev- 
eral observed anomalies in the pseudogap phase of the 
cupratea^i^i^iSiLSiSi. Attempts to directly observe it in 
neutron scattering are encouragingiS. 

One of the arguments against the DDW theory is that 
it is inconsistent with in-plane optical conductivity mea- 
surements. The assumption has been that as the tem- 
perature is reduced below the pseudogap transition tem- 
perature T*, and a (partial) gap opens, some of the low- 
energy excitations are lost and this should cause a re- 
duction in the low-frequency optical conductivity, which 
is compensated by an increase in the spectral weight 
at higher frequencies. Such behavior is seen in some 
charge density-wave systemsii, but is not observed in the 
cupratesi2iiii4, which might be taken as an indication 
that the pseudogap is associated with superconducting 
fluctuations, which would cause a downward motion of 
spectral weight (which would eventually coalesce into a 
zero-frequency (5- function in the superconducting state). 

These experiments find no changes in the in-plane 
spectral weight within their error barsiSiiii^: spectral 
weight lost at the low-energy end of the measured fre- 
quency range appears (perhaps unexpectedly) at very low 
frequencies, and hence accompanies a narrowing of the 
low-energy optical conductivity peak. This phenomenon 
was interpreted in ref. ^3 ^ the lack of any in-plane 
optical evidence for the pseudogap. 

In this paper we show how the DDW theory of the 
pseudogap can be consistent with this observed effect. 
The key to this phenomenon, we believe, (as conjectured 
in ll2|) is the strong anisotropy in the quasiparticle spec- 



trum and scattering rate in the cuprates. ARPES exper- 
iments clearly show that quasiparticles in the antinodal 
region of the Brillouin zone {{kx, ky) ~ (0, ±7r), (±7r, 0)), 
where the gap opens in the pseudogap phase, scatter 
much more strongly than they do in the nodal region 
{{kx,ky) ~ (±7r, ±7r)), where the gap is zero. (While 
the transport lifetime is not the same as the quasiparti- 
cle lifetime observed in ARPES, we assume that it has 
a similar anisotropy, as in "hot spot" and "cold spot" 
theoriesiSiii. ) The lost carriers, therefore, give a rela- 
tively small contribution to the conductivity. The loss of 
these carriers can be cancelled by a general reduction of 
all the scattering rates, as the temperature is reduced. 
As pointed out in ref. Ila . this may also explain the ob- 
served behavior in several classes of spin-density-wave 
(SDW) and CDW systems. 

These ideas are borne out by a model calculation with 
the mean-field DDW Hamiltonian and an ansatz for the 
quasiparticle lifetime. We show that when the lifetime 
has an anisotropic form motivated by ARPES measure- 
ments, the conductivity spectral weight shifts downward, 
as in the cuprates. However, for an isotropic lifetime, the 
spectral weight shifts upward. 



II. MEAN-FIELD HAMILTONIAN 

The mean-field Hamiltonian for the DDW state is 

-ffDDW = X![(^k-M)4^Ck^ + («Akcj^^Ck+Qff+ h.c.)] (1) 
k,(T 

where Ck is the annihilation operator of an electron in a 
state with momentum k and spin a, fj, is the chemical 
potential, Ak = (cj^^rCk+Qcr) = Ao(cosfc:c - cosfcj^)/2 is 
the DDW order parameter and Q = (tt, tt) is the DDW 
ordering wave vector, with lattice spacing set to unity. 
The tight-binding band structure is given by ek — eik + 
e2k, with 

eik = — 2t(cosfc2; -|-cosfcj^) , £2)1 ~ 4:t' cos kx cos ky. (2) 



2 




FIG. 1: DDW quasiparticle energy bands along the line ky — 
in the Brillouin zone for t = 0.3, t' = 0.09, Ao = 0.2 and 
fj, = -0.3 eV. 



where t and t' are nearest-neighbor and next-neighbor 
hopping parameters. 

Introducing the two-component field operator xl^a ~ 



kcr '■^ 

written as 



c^^ ^*"^k+Qcr I mean field Hamiltonian can be 



HdDW — ^ Xkcr-SkXko- 



with 



Ek - M ^k 
^k Ck+Q — ^■ 



Bk = (e2k - m) + eikcr^ + Akcr^ 



(3) 



(4) 



(5) 



where cr' are the PauH matrices and the sum is over 
half the original Brillouin zone (reduced Brillouin zone 
- RBZ), i.e. \k^\ + \ky\ < tt, . 

Diagonalizing the Hamiltonian will then give the DDW 
quasiparticle energy bands 



Eia = (e2k - /i) ± Ve?k + ^k 



(6) 



as depicted in figure ^ 



the only relevant interactions are the electron-electron 
interactions which generate the DDW coupling. All 
other interactions including quasiparticle-quasiparticle 
and quasiparticle-impurity interactions will be taken into 
account later by assuming a non-zero self energy). 

The 2x2 Nambu Green function is defined by 



Go(k,t)-(rxk(<)x|,(o)) 

^ / (Tck(t)4(0)) »(Tck(t)c^Q(0)) 

1^ -Z(rck+Q(i)4(0)) (TCk+Q(t)4+Q(0)) ^ 

The Fourier transform of the Green function matrix is 
then 



(7) 



(w -I- iS) - (e2k - m) - eikcr^ - AkCT^ 

_ - (£2k - + eik^^ + AkCr^ 

.2 , A2 
-Ik ^ ^k 



(8) 



We now consider the effects of impurity scattering and 
'residual' electron-electron interactions. Here, we will 
capture their combined effect by introducing a non-zero 
single-particle self-energy I](k, uj) = Si(k, a;)-|-iE2(k, lu), 
where the real and imaginary parts give the energy renor- 
malization and quasiparticle lifetime, respectively. Ne- 
glecting the shift in the excitation energy due to the real 
part Si(k, cj), the self-energy can be written in terms of 
the quasiparticle lifetime as 



E(k,w) 



T(k,C.) 



(9) 



where l/r(k, w) is the quasiparicle scattering rate (in- 
verse lifetime) . Including the self-energy, the full Green's 
function G(k, uj) will then be given according to Dyson's 
equation by 



III. GREEN FUNCTION 



G(k,a;) = (Go(k, 



^ik,uj)y 



(10) 



The "non- interacting" Nambu Green's function can 
then be obtained by inverting the matrix B^ . (For now. 



The spectral function A(k, uj) can then be calculated 
by taking the imaginary part of the Green function 



J 



A{k,uj) = Im G(k,w) 

TT TTT 



1 (lo - e2k + ^^)^ + (e?k + Ak) + + 2(w - esk + ^J-)ielk(T^ + Ak^^) 



[{lu - Ei^)iuj - E2k) - (l/r)2]2 + (2{uj - e2k + m)/^)^ 
I 



(11) 



In Figure[21 the diagonal elements of the Nambu spec- tral function matrix A{uj, k^, ky = 0) are plotted against 
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with 



n,.(T) = (T,j,(r)j,(0)) (14) 
The current operator for the DDW quasiparticles can 
be obtained by minimaUy-coupling the mean-field Hamil- 
tonian |^ to the electromagnetic field, A, and differen- 
tiating once with respect to A: 



j = J2 Vf2(k) (xlxk) +VFi(k) (^xWxk 



RBZ 



VA(k) (xi-Cr^Xk 



(15) 



where Vi^i(k) = Vkei(k) ,VF2(k) = Vke2(k) and 
VA(k) = VkA(k). 

Using this form of the current operator, the current- 
current correlation function can be written in terms of 
the elements of the imaginary-time Nambu Green's func- 
tion Qij (the non-interacting Green's function Gy is now 
promoted to the interacting one Qij , to take the effect of 
the residual interactions into account). We will have 



FIG. 2: Elements of the DDW spectral function matrix for 
1/r = 0.005 eV: {a.)Aii{uj,kx,ky = 0) and (h)A22{i^,kx,ky = 
0). Other parameters are the same as in figure Q 



kx and uj. For demonstration purposes, the scattering 
rate here is assumed to be a constant 1/r = 0.02 eV (in- 
dependent of momentum and frequency). The diagonal 
entries An and A22 have peaks centered mostly at ener- 
gies corresponding to the upper and lower energy bands 
{El and E2) respectively. 



{Trj{r)jm = - ^ {[^F2mhr{g{-r)g{r)) 



RBZ 

+ [vpi{k)]Ma'g{-T)a'g{T)) 

+ [vA(k)]2tr {<j'g{-Tyg{T))) 

+ (vFi.VF2)tr ((fT3g(-r) + gi-T)a^)g{T))) 
+ (vf.2.VA)tr {ia'g{-T) + g{-T)a')g{T))) 
- (vi.i.VA)tr {{a^g{-T)a'' + a''g{-T)a^)giT))) 



(16) 



IV. OPTICAL CONDUCTIVITY 

Having found the spectral function A{Lu,k), the real 
part of the AC conductivity can now be calculated by 
using the Kubo formula 

Re a^xi^^) = —'^m Uxxiii^n + iS) (12) 

OJ 

where Il{iujn) is the Fourier tansform of the current- 
current correlation function in Matsubara formalism 

n,,(zc^„)= / dre'^-^nxxir) (13) 
Jo 



In this equation, we have ignored vertex corrections. 
These are important when the scattering rate is strongly 
angle-dependent, as they distinguish the transport and 
quasiparticle lifetimes (for instance, through a (1 — cos 61) 
factor) and also distinguishing umklapp scattering from 
momentum-conserving scattering. In what follows, we 
will assume that the replacement r — > rtr is made {i.e. 
ignore vertex corrections). Referencei^ has considered 
the vertex correction for DDW conductivity. 

Writing g in terms of the spectral function A{]i,uj), 
evaluating the Matsubara sum, and doing the analytic 
continuation {iujn — * co + iS), we find the optical conduc- 
tivity to bei^ 
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(^) _ i ^ /" ^\nF{e) - nF{e + uj)] 



e + uj) + 2Ai2{k, e)Ai2{k, e + lu) + A22{k, e)A22{k, e + uj) 
+ [vFi(k)]^ An{k, e)Aii(k, e + uj) - 2Ai2{k, e)Ai2{k, e + uj) + A22{k, e)A22{k, e + uj) 
ik)f\A22{k,e)Anik, e + uj) + 2Ai2{k, e)Ai2{k, e + uj) + Aii(k, e)A22{k, e + uj) 
+ 2vFi(k).VF2(k) Aii{k, e)Aii{k, e + uj) - A22{k, e)A22{k, s + uj) 
+ 2vF2 (k) .VA (k) Ai2 {k, e) {Aii{k,e + u) + A22{k,e + u)) + {An {k, e) + A22 {k, e))Ai2 {k,e + uj) 
+ 2vFi(k).VA(k) \Ai2{k,e){Aii{k,e + uj)- A22{k,e + uj)) + {Aii{k,e) - A22{k,e))Ai2{k,e + u) 



(17) 



T =0.03; A,=0.25; 
y =0.02; 




FIG. 3: Real part of the optical conductivity for a constant 
scattering rate 7 = 70 = 0.02, Ao = 0.25 and T* = 0.030. 



where nF{e) is the Fermi distribution function. For 
demonstration purposes, in fig. |31 the real part of the 
optical conductivity has been plotted against uj for two 
different temperatures (one above T* , depicted with a 
solid line, one below T*, depicted with a dashed-dotted 
line), assuming that the quasiparticle lifetime is a con- 
stant (temperature and momentum independent) and the 
gap is unrealistically big (Wq = 0.25 eV). As expected an 
upward shift of the SW occurs when the gap opens. Sim- 
ilar calculations have also been done in2&. 



ticularly in the anti-nodal regions, where even lowest- 
order perturbation theory around the DDW mean-field 
Hamiltonian^- predicts short lifetimes which may indicate 
a breakdown of quasiparticles. However, we will compute 
the conductivity in the quasiparticle approximation to 
show that, even at this level, an upward shift of spectral 
weight is not expected. If there are no quasiparticles at 
the anti-nodes, then the situation may be even better. 

In order to proceed with this strategy, we need one final 
ingredient, an ansatz for the quasiparticle scattering rate 
1/T(k, u;;T) = 7(k, ti;;T), as a function of momentum, 
frequency, and temperature in the underdoped cuprates. 
A number of angle-resolved photoemission experiments 
have measured the inverse lifetime (imaginary part of the 
self energy), as a function of these different parameters. 
In these experiments, the width of the quasiparticle peaks 
in the energy distribution curves (EDC's) or momentum 
distribution curve (MDC's) are measured as functions 
of momentum, energy and temperatureSLSSiSiSi. While 
the transport lifetimes are not necessarily identical to the 
quasiparticle lifetimes (or, equivalently, the vertex correc- 
tions are not necessarily small) , we expect that they will 
have a similar anisotropic behavior. For the purposes 
of this model calculation, we will simply take them to 
be the same. We emphasize that the lifetimes which we 
adopt below are for illustrative purposes since our main 
goal is to show that an upward movement of spectral 
weight is not a necessary concomitant of the emergence 
of DDW order at finite temperature. We are not making 
any claims here about the correctness of these lifetimes. 

We take the imaginary part of the self-energy (quasi- 
particle scattering rate) to have the form 



V. QUASIPARTICLE LIFETIME 

Applying the formulas of the preceding section to the 
underdoped cuprates presupposes that the quasiparticle 
picture makes sense there. This is questionable, par- 



I]2(k,c^;T) = S2(c^;T)-t-r(k). 



(18) 



where Yi2{^]T) is temperature and frequency dependent 
with no momentum dependence and r(k) is strongly mo- 
mentum dependent. 

Such a form has been motivated by marginal Fermi 
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liquid phenomenolog}i2SiS& (MFL) . In this way of analyz- 
ing the data, it is assumed that the quasiparticle lifetime 
which comes from electron-electron scattering is indepen- 
dent of temperature and linear in energy for small tem- 
peratures and linear in temperature, independent of the 
binding energy for higher temperatures: 



^MFL 



{uj;T) = AMax(|w|,T). 



(19) 



In this ansatz, all of the angular dependence comes from 
elastic electron-electron scattering. This is a convenient 
form, but we will show that our results hold even for 
some others. 

For instance, we repeat our calculations with the stan- 
dard Fermi liquid (FL) quasiparticle lifetime: 



(20) 



again assuming that the angular dependence comes from 
r(k). In order to show that these particular forms of uj- 
dependence are not playing a big role in shifts of SW, we 
have also tried w-independent forms of these scattering 



T— Linear 



(w; T) = AT and T,^ (w; T) = AT^. 



rates: E 

While this form of the quasiparticle lifetime does not 
give the correct DC conductivity, it is still useful as 
a check because our goal is to emphasize the role of 
anisotropy and show that it can lead to a downward shift 
of spectral weight, regardless of its detailed frequency and 
temperature dependence. 

The quasiparticle lifetime is strongly anisotropic^: ex- 
citations in the antinodal region arc more strongly scat- 
tered than the ones in the nodal region by up to a factor 
of 5. Hence, the assumed form (|18|) necessitates that 
r(k) be a strongly anisotropic function of k. We take 
the form: 



rAniso(k) = 7o(l + (C0sfc:r " COS fcj^ 



(21) 



where 70 is the scattering strength in the nodal region. 

In order to check the role of anisotropy in the eventual 
shifting of the SW, we also check the case in which r(k) 
is (unrealistically) isotropic: 



riso(k) = 70. 



(22) 



a: 



(a) 



>. = 0.1, y,=0.02 

1=0.030" 

T=0.001 





FIG. 4: Real part of the optical conductivity for (a) isotropic 
scattering rate with 70 = 0.020 and (b) anisotropic scattering 
rate with 70 — 0.008. Quasiparticle lifetime is taken to be 
linear in temperature and independent of u. 



VI. RESULTS 



In Figures 01 to [7| we have plotted the real part of 
the optical conductivity vs w, when the momentum- 
independent part of the lifetime is given by the four dif- 
ferent forms listed above. In each figure the isotropic 
case is compared with the anisotropic one with the same 
set of parameters (listed in the figure captions) except 
for 7o , which is smaller in the anisotropic case because of 
its extra (larger than 1) prefactor. In all cases we have 
t = 0.3, t' = 0.09, /i = -0.3 and T* = 0.03. 

It is clear that in figures 4a, 5a, 6a, 7a, in which the 
scattering rate is isotropic, there is an upward movement 



of spectral weight (though it is small in some cases). 
However, in figures 4b, 5b, 6b, 7b, in which the scatter- 
ing rate is anisotropic, there is a clear downward move- 
ment of spectral weight. (Broadening the Drude peak 
has also been seen in22i andi^, in which scattering has 
been isotropic.) 
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(a) 

■f = y + X Max ((0 , T) 



x = o.^, Y=o.o2 



T=0.030 
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(a) 



1.0, Y,=0 01 

T=0.030 

— - T=0.001 



0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 




FIG. 5: Real part of the optical conductivity for (a) isotropic 
scattering rate with 70 = 0.020 and (b) anisotropic scattering 
rate with 70 — 0.008. Quasiparticle lifetime's temperature 
and frequency dependence is given by the marginal Fermi liq- 
uid theory. 



FIG. 6: Real part of the optical conductivity for (a) isotropic 
scattering rate with 70 = 0.010 and (b) anisotropic scattering 
rate with 70 — 0.002. Quasiparticle lifetime is taken be to 
independent of lo have a temperature dependence. 



VII. CONCLUSIONS 

Recent data by Santander-Syro et a'M., has questioned 
the previous belief that the opening of the pseudogap in 
the underdoped cuprates can be detected by in-plane op- 
tical conductivity measurements. The unexpected result 
that as the temperature is reduced (and although a gap 
opens in the antinodal region of the Brillouin zone), the 
spectral weight is still transferred to lower frequencies, 



was interpreted as lack of any pseudogap signature in 
the optical data. In this paper we showed that this ef- 
fect is consistent with the DDW theory of the pseudogap 
state of the underdoped cuprates. 

In the four sets of graphs of the previous section it 
can clearly be seen that the key factor in deciding which 
way the SW is transferred is isotropy or anisotropy of 
the scattering rate. It is shown that regardless of the 
form of the temperature and frequency dependence of 
the quasiparticle lifetime {e.g. Fermi liquid or non- Fermi 
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y = Y { 1 + (Cos kx - Cos kyf )+ Max (to" , T") 



= 1.0, 
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T=0.001 



liquid) the SW is shifted upward for the isotropic scatter- 
ing rate while it is shifted downward for the anisotropic 
case. (We have not tried to find a form for the quasiparti- 
cle lifetime which correctly reproduces all of the transport 
data, but have focussed on the role of anisotropy and have 
shown that similar results are obtained for several differ- 
ent forms of frequency and temperature dependence.) 



We believe the explanation to be as follows: the pseu- 
dogap opens in the antinodal region where the carriers 
are more strongly scattered than the ones in the nodal 
region, where there is no gap. Therefore, we have lost 
those excitations which already gave relatively little con- 
tribution to the transport properties of the normal state. 
Furthermore, by reducing the temperature, the scatter- 
ing rate of all the excitations is reduced. Our results 
clearly show that in the (more realistic) anisotropic case, 
the effect of lost excitations due to the gap opening can 
be more than canceled by a temperature-dependent re- 
duction of the scattering rate for the rest of the excita- 
tions and hence a downward shift of the optical spectral 
weight. It is obvious that anisotropy is the key here, 
since for similar parameters for isotropic scattering rate, 
an upward transfer is observed. We note that similar 
considerations should apply to the case of 2H-TaSe2^. 



FIG. 7: Real part of the optical conductivity for (a) isotropic 
scattering rate with 70 = 0.020 and (b) anisotropic scattering 
rate with 70 = 0.008. Quasiparticle lifetime's temperature 
and frequency dependence is given by the Fermi liquid theory. 
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